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Abstract 

This paper is a slightly modified version of the introductory part of a PhD thesis, also 
containing the articles hep-ph/0303019, hep-ph/0510375 and hep-ph/0512177 We provide 
a short history of the research of electroweak thermodynamics and a brief introduction to 
the theory as well as to the necessary theoretical tools needed to work at finite tempera- 
tures. We then review computations regarding the pressure of electroweak matter at high 
temperatures (the full expression of the perturbative expansion of the pressure is given in 
the appendix) and the electroweak phase diagram at finite chemical potentials. Finally, 
we compare electroweak and QCD thermodynamics. 
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Chapter 1 

Introduction 



Understanding the properties of matter is one of the main goals of modern physics. Under 
normal conditions, matter is composed of atoms, molecules and free electrons and interactions 
between them can be described by (quantum) electrodynamics. However, when the temper- 
ature and density of matter are increased, these basic building blocks of matter as we know 
them begin to break apart and matter eventually becomes a collection of elementary particles. 
Then we need to resort to theories of particle physics to describe its properties. 

The elementary particles as well as the strong, weak and electromagnetic interactions 
between them are to high accuracy described by the standard model of particle physics. It 
contains a rich collection of interesting features at extreme conditions, related to, for example, 
different phase transitions. At present, only the properties of strongly interacting matter, 
QCD plasma and hadron gas, are within the reach of experimental studies. Knowledge about 
the behavior of electroweak matter is, on the other hand, purely theoretical at the moment. 
However, it is important that we understand properties of weakly interacting matter as well, 
since it is suggested by modern cosmology that in the very early universe temperature was so 
high that such matter existed. 

In this thesis we will study some aspects of electroweak matter. The thesis consists of this 
introductory part and of three research publications [1-3] . The first paper [1] considers the 
electroweak phase diagram. In it, high temperature dimensional reduction of the standard 
model was formulated in the presence of non-zero chemical potentials for baryon and lepton 
numbers and the computation was then applied to calculate the phase diagram. Special 
attention was given to the location of the endpoint of the first order phase transition line. 
The papers [2, 3] are concerned with calculating the pressure of electroweak matter at finite 
temperature and zero chemical potentials. In Ref. [2] the pressure was calculated at high 
temperatures to three loops, or to order g 5 in the coupling constants and its properties, such 
as scale dependence, were analyzed. Comparison with QCD pressure was also performed. In 
Ref. [3] the previous calculation was extended to lower temperatures (temperatures around 
the electroweak crossover) by reorganizing the effective field theories used to calculate the 
pressure. 

This introductory part is organized as follows. We will first review some history of elec- 
troweak thermodynamics and then in chapter 2 present the basic structure of the electroweak 
theory. In chapter 3 we consider field theories at finite temperatures and chemical potentials 
in general. The actual calculations concerning the pressure of the standard model (chapter 
4) and electroweak phase diagram (chapter 5) are discussed next. In the final chapter we 
consider the similarities and differences between the QCD and electroweak thermodynamics 
relevant to the thesis. Expansion of the pressure of the standard model to three loop order is 
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given in detail in the appendix A. 



1.1 History 

Soon after the electroweak model was constructed, interest about its consequences when ap- 
plied to thermal systems arose. Based on the close analogy between the bosonic sector of 
the electroweak model and the Ginzburg-Landau theory of superconductivity, Kirzhnits and 
Linde proposed that the symmetry that is spontaneously broken at low temperatures would be 
restored at high temperatures [4,5]. The reason is that, unlike at low temperatures where the 
equilibrium state of the system is such that energy is minimized, at high temperatures thermal 
equilibrium is achieved when entropy of the system is maximized, 1 i.e., when the symmetry 
of the theory is restored (no order in the system). They estimated that the restoration of 
symmetry would happen in a second order phase transition at temperatures of the order of 
T ~ G^ 1 ^ 2 ~ 10 2 GeV, where is the Fermi coupling constant. A more systematic approach 
using finite temperature effective potentials was then developed by Dolan and Jackiw [6] and 
by Weinberg [7]. Those studies confirmed the validity of earlier, more heuristic, arguments 
and a more quantitative estimate for the critical temperature was derived for a number of 
different theories. Electroweak phase transition thus gained a status as a basic ingredient in 
cosmology. The possibility that the transition would be of first order was shortly after taken 
into account by Kirzhnits and Linde [8] who calculated leading order corrections to the pre- 
viously evaluated effective potentials. Similar phase transitions in grand unified theories lead 
to the development of inflationary models in cosmology [9]. 

Need for a more detailed understanding of the electroweak symmetry restoration became 
relevant when it was realized that the electroweak phase transition might have a crucial role 
in understanding the baryon number asymmetry in the universe [10, 11]. Electroweak vacuum 
has a non-trivial topological structure and transitions between topologically distinct vacua 
violate conservation of baryon number which could lead to erasure of any baryon number 
asymmetry in the universe. On the other hand, the same transitions could also provide a 
mechanism for producing the baryon number asymmetry, depending on the properties of the 
electroweak phase transition [12]. If the transition is strongly first order, a baryon number 
asymmetry could be generated during the transition. The prospect of understanding the roots 
of the baryon number asymmetry in the universe lead to a renewed interest in quantitative 
description of the electroweak phase transition. 

Building on the work laid out by Dolan, Jackiw and Weinberg, the electroweak phase 
transition was first studied in detail with one-loop effective potential calculations [13-15]. At 
high temperatures and assuming that the Higgs boson is sufficiently light, the strictly one- loop 
expansion of the effective potential is given by [13] 

V(<p, T) = D(T 2 - T V - ET^ + , (1.1) 

where (p is the expectation value of the Higgs field and the different coefficients can be com- 
1 Wc arc considering canonical ensembles. 



2 




Figure 1.1: A schematic plot of the evolution of the scalar potential as temperature is decreased 
below the critical temperature. 



puted to be 



D = (2m 2 w + m 2 z + 2m 2 t ) , 

E = ^(2rn* w + rn%), 

A, = A-^( 2 ^ln^ + m |ln-^-4^1n-^), 

B = 64^ (2m4vK + m4z " 4m ^' 

3 3 
In c B = - + 21n4vr-27, In c F = - + 21nvr - 2j. (1.2) 

Here all the masses are measured at zero temperature and ipo = 246 GeV is the value of the 
scalar condensate at T = 0. 

At high temperatures the minimum of this potential is achieved when the expectation 
value of the scalar field vanishes, cp = 0, and thus the symmetry is exact. However, as the 
temperature is decreased, another local minimum appears, which becomes the global minimum 
at some critical temperature T c and it becomes favorable for the system to reside there (see 
Fig. II. 1|) . Thus, the symmetry of the theory gets spontaneously broken. Using the expression 
for the effective potential above, it is possible to calculate many of the essential features of 
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this phase transition. For example, the critical temperature is easily evaluated to satisfy 



Another important quantity is the value of the expectation value of the scalar condensate in the 
broken symmetry phase at the moment of the phase transition, ip c = 2ET c /Xt c . Baryogenesis 
can be explained within electroweak physics only if the ratio y? c /T c is large enough. The 
baryon number violating processes must be cut-off after the transition or any baryon number 
asymmetry that was generated in the transition will be washed out. Sufficient criterion for 
this is that ip c /T c > 1 is fulfilled [12]. 

Although the one-loop approximation can be used to calculate the characteristics of the 
phase transition, it is not guaranteed to be a reliable method. Indeed, perturbative calculations 
in gauge field theories are known to suffer from infrared problems at high temperatures [16,17]. 
Collective phenomena such as screening of electric fields (Debye screening) , not contained in 
the strict loop expansion, take place. The Debye screening can be taken into account by 
introducing a thermal mass for the static Coulomb fields, ~ gT. This corresponds to 
resumming an infinite class of ring diagrams to the effective potential and thus yields an 
improved one-loop expansion of the effective potential [14,15]. The essential difference to the 
strict one-loop expansion of the potential is that the value of the coefficient E of the ~ y? 3 
term, responsible for the first order nature of the transition, will be reduced by a factor of 2/3 
and thus the transition becomes weaker. Especially, the result suggests that the transition will 
not be strong enough to explain the baryon number asymmetry if the Higgs mass is greater 
than m H > 45 GeV [15]. 

The analysis has been further improved by calculating the two-loop corrections to the 
effective potential, to order g 4 ,X by Arnold and Espinosa [18] and to order g 4 ,X 2 by Fodor 
and Hebecker [19]. Although the two-loop calculation provides (at least formally) a more 
precise description of the phase transition, the qualitative features remain the same as in the 
one-loop calculation, predicting a first order phase transition. 

With a quantitative understanding provided by the effective potentials, it is possible to 
study the dynamics of the electroweak phase transition in more detail [13,15,20-22]. A first 
order phase transition proceeds by nucleation of true equilibrium state bubbles into the system 
which subsequently grow until the whole universe resides in the true equilibrium state. Such a 
nucleation process leads to local departure from thermal equilibrium and can leave observable 
traces in the universe. 

The studies described above, based on calculating the properties of the phase transition 
with perturbatively computed effective potentials, though convenient, are ultimately unreliable 
due to infrared divergences that plague perturbative gauge field theories [16,17]. Even though 
the infrared problems are not manifest when calculating in the broken symmetry phase (since 
particles obtain masses via Higgs mechanism which serves to regulate the infrared divergences) , 
to obtain information about the phase transition one must compare the calculations in the 
broken symmetry phase to those in the symmetric phase where the infrared problems reappear. 
To obtain reliable results, one must therefore resort to non-perturbative methods. 

The non-perturbative analysis is convenient to perform by identifying the modes that are 
responsible for the infrared divergences in perturbation theory. These modes can be isolated 
to a series of effective theories [23-28] and we can combine the use of perturbative calculations 
(applied to infrared safe modes) and numerical computations (applied to modes that are 
infrared divergent in perturbation theory). Carrying out such calculations, the phase diagram 



To 



(1.3) 




4 



The Standard Model 
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Figure 1.2: The phase diagram of the electroweak theory. 



of the electroweak theory was solved by a number of authors [29-33]. What is observed (see 
the phase diagram in Fig. 11.21 given in [34]) is that, although for small Higgs masses there 
is a first order phase transition in the electroweak theory, as the Higgs mass becomes larger, 
the phase transition weakens. The first order phase transition line has a 2nd order endpoint 
at rriH ~ 72 GeV and for larger Higgs masses no phase transition is observed, just a smooth 
crossover. Furthermore, the endpoint has been observed to be of 3d Ising universality class [35]. 
Since there is no first order phase transition in the (minimal) electroweak theory for physical 
values of the Higgs mass, electroweak baryogenesis (within minimal standard model) has been 
ruled out. 

In addition to considering just the effects of high temperatures, properties of electroweak 
matter have also been studied when density of matter is high. It was noted by Linde [36], 
using again analogy to superconductivity, that increasing fermion number density would lead 
to an increase in the symmetry breaking, i.e., the expectation value of the scalar condensate 
would increase (conversely, non-zero fermion current density tends to restore the symmetry). 
These considerations were brought to a firmer ground by computing the one-loop effective 
potential with finite chemical potentials associated with conserved fermionic numbers taken 
into account [37-41]. Similarities between spontaneous symmetry breaking and Bose-Einstein 
condensation by considering finite bosonic chemical potentials (related to conserved gauge 
charges) have also been studied [42,43]. 

Complete thermodynamic description of the standard model is given when, in addition to 
temperature and chemical potentials for the conserved particle numbers, also the effects of an 
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external U(l) magnetic field are taken into account. It is much speculated that strong magnetic 
fields might have been present in the very early universe, seeding the galactic and inter-galactic 
magnetic fields that we observe today [44]. The effect that such magnetic fields could have 
on the electroweak phase transition has been studied [45] and it was observed that they can 
strengthen the transition somewhat and lower the transition temperature, but not enough 
to save electroweak baryogenesis. It was also noted that no exotic phase, analogous to the 
Abrikosov vortex lattice of type-II superconductors and implied by perturbative calculations 
[46], could be found. 
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Chapter 2 



Basic structure of the electroweak 
theory 

The standard model of particle physics is a quantum gauge field theory based on local SU(3)c x 
SU(2)i x U(l)y symmetry. It describes dynamics between leptons (SU(3)c singlets) and 
quarks (belonging to the fundamental representation of SU(3)c). These are divided into 3 
families, each composed of 2 leptons and 2 quarks. The families are identical to each other 
in every other respect apart from the Yukawa couplings giving masses to the particles. This 
theory has been confirmed experimentally with an extraordinary precision and it forms the 
cornerstone of modern particle physics. 

In this thesis we will concentrate on the electroweak sector of the standard model. The 
electroweak theory was first formulated by Weinberg, Glashow and Salam [47-49] and it 
replaced Fermi's theory 1 of /3-decay as a description of weak interactions. It predicted neutral 
current processes, unknown at the time, and the discovery of those [50] and of the weak gauge 
bosons [51] established the status of the theory. 

In this chapter we will briefly review the basic setting of the electroweak theory and its 
properties relevant to the thesis. The theory is defined by the Lagrangian 

— -^bos. ~i~ -^ferni. ~i~ ^QCD ^top mass; 

where the different parts are defined in the following. The gauge fixing and ghost terms in 
the Lagrangian are not explicitly written down. 

2.1 Higgs mechanism 

A distinct feature of electroweak interactions is that the and Z° bosons mediating the 
interactions are massive. Since gauge invariance protects gauge bosons from acquiring masses, 
it would then seem impossible to describe weak interactions in terms of a gauge field theory. 
However, even though the theory has a symmetry, it is not necessary that the ground state of 
the theory has the same symmetry, that is, the symmetry may be spontaneously broken. This 
is a sufficient requirement for producing masses for gauge bosons. In the standard model, this 
is accomplished by introducing a scalar field, called the Higgs scalar, into the theory. The 
scalar belongs to the fundamental representation of the SU(2) symmetry group and acquires 
a non-zero vacuum expectation value that serves to spontaneously break the symmetry. 

1 \t was already known before the electroweak theory was formulated that the Fermi coupling could be 
understood as an exchange of massive, charged vector bosons between fermions. 
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Consider the bosonic sector of the electroweak theory. The Lagrangian describing the 
interactions is 



£ bos . = -\G%G^ a - jF^F^ + D^D^ + v 2 &<& - A (V$) 2 , (2.2) 

where <E> is the scalar doublet in the fundamental representation of SU(2), D^Q = — 
igA^T 11 ^ /2 — Yig'B^Q with r a being the Pauli spin matrices and Y = 1/2 is the hypercharge 
of the scalar doublet. The field strength tensors are given by G" = d^A^ — d u A a ^ + ge abc 'A b A c v 
and F^y = d^By — d u B^ where A°^ and B^ are the gauge bosons of weak and hypercharge 
interactions, respectively, and e abc are the generators of the adjoint representation of SU(2). 
The symmetric vacuum <3? = is unstable due to the sign of the mass term for the scalar 
doublet (see Fig. I2.1j) . At tree level, a stable solution to the equations of motion satisfies 
instead $o$o = ^ 2 /(2A) and is commonly chosen to be 

where v 2 = v 2 /A. We should consider fluctuations around this state. The vacuum is no longer 
invariant under the full gauge group, thus spontaneously breaking the symmetry. A residual 
U(l) symmetry, generated by Q = r 3 /2 + Y, remains exact. The mass eigenstates and the 
masses can be readily worked out 



W, 



4 = mW-^). m l = W+9 2 )v\ (2.4) 

V 9 +9 4 

A * = HT7-^[ (9'Al + gBj, m 2 = 0. 

V 9 + 9 

The photon A„, corresponding to the residual U(l) symmetry, remains massless as it should. 
The electric coupling e is related to the coupling constants g and g' by e = gg' / \J g 2 + g' 2 . 

Although it is experimentally verified that the electroweak gauge bosons are massive, full 
confirmation of the Higgs mechanism is still lacking since the Higgs boson has not been found. 
Today, the direct experimental lower limit for the Higgs mass is mu > 114 GeV [52]. In spite 
of this, we will assume in this thesis that the Higgs mechanism of the minimal standard model 
is valid. 



2.2 Fermionic sector of the theory 

Unique feature of electroweak interactions is that the different chiral projections of fermions 
transform differently under the gauge group. The left handed components form doublets under 
SU(2) transformations while the right handed components are SU(2) singlets. This leads to 
parity violation in electroweak processes. The fermions are classified into families as 



leptons 
quarks 
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Figure 2.1: A schematic plot of the Higgs potential. 



where the subscripts L/R refer to left /right handed components, defined by 

ipL = -(1-75)^, tpR = -(l + 75)-0- ( 2 - 5 ) 

Note that there are no right handed neutrinos in the minimal standard model. Denoting by 
II and e R the left handed doublets and right handed singlets of leptons, respectively, and by 
qL, ur and da the left handed quark doublets and right handed u and d -type quark singlets, 
respectively, the coupling of fermions to the gauge fields is governed by 

£fcrm. = y^&Pj), ip S {lL,e R ,q L ,u R ,d R }, 

Dp = d^-IigA^-Yig'B^, (2.6) 

where / and Y are the weak isospin and the hypercharge of the relevant fermion, respectively. 
The values of those are given in table l2~Tl 

Due to the chiral transformation properties of fermions, introducing masses for them using 
standard mass terms mipip is not possible since such terms would violate gauge invariance. 
Instead, masses can be generated for them by coupling them to the Higgs field through Yukawa 
couplings. Since all the fermions except for the top quark are so light that their masses can 
be neglected for our purposes, we only need to take into account the top Yukawa coupling 

Aopmass = igyq^^tR + h- c - (2.7) 
In the symmetry breaking the mass nit = 9yv/V2 is generated for the top quark. 
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QL 


1/2 


1/6 


2/3,-1/3 


UR 





2/3 


2/3 


d-R 





-1/3 


-1/3 



Table 2.1: Values of the various charges for fermions. Here T 3 stands for the third component 
of weak isospin. For doublets, the electric charge is given for each member of the doublet. 



Since the theory contains quarks, we have to take the strong interactions between them 
into account. These are governed by QCD, 

£qcd = -\w^W^ a -ig s ^C-T a q, (2.8) 

quarks 

where q are the quarks, W« v = d^C* - d v C« + g s f abc C b C^ with C£ being the gluons, T a and 
jabc g enera tors of the fundamental and adjoint representations of SU(3)c, respectively 
and g s is the strong coupling constant. Note that the free propagation of quarks, qicjtq, is 
already taken into account in £f erm . 

2.2.1 Conservation of fermion numbers 

At classical level, the baryon number and the three lepton numbers are conserved indepen- 
dently. 2 Instead, in the quantum theory, due to the chiral couplings of fermions to the weak 
gauge bosons, there is an anomaly [53-55] and the conservation of baryon and lepton numbers 
is broken, 3 

^if = ^wTrG^G'* 7 , (2.9) 

where the jjf correspond to the baryon and the lepton number currents. The baryon and lepton 
number violations are, however, related so that (A£>)/3 = AL e = AL^ = AL T . Consequently, 
the three linear combinations 

Xi = -B-Li, i = e,v,r (2.10) 

are exactly conserved, while the remaining combination, B + L, is not. Here B stands for 
the baryon number, Lj for the lepton numbers in each family and L = £^ Lj. From now on, 
we will refer to the non-conserved B + L number as "baryon number" and to the conserved 
B/3 — Li numbers as "lepton numbers". 

From equation (|2.9|) we see that in order for B + L to change, the gauge fields must evolve 
so that (for a review on processes governing the baryon/lepton number non-conservation, 

2 The baryon numbers for each family are not separately conserved since for quarks, the weak eigenstates are 
not equal to the mass eigenstates. If we considered a model where neutrinos are massive, then similar mixing 
would be possible also in the leptonic sector and thus there would just be one conserved lepton number. 

3 Even though the baryon and lepton number currents are not conserved, the charge currents are and hence 
there is no anomalous breaking of gauge invariance, which would spoil renormalizability of the theory. 



10 



see [56]) 



A(B + L) = ^Lj\ t Jd 3 xe^ pa TrG^G^^0. (2.11) 

In general, this requires very strong fields and thus such processes are not present under 
normal conditions. The non-conservation of B + L is, however, relevant when one considers 
the vacuum structure of the theory. The vacuum has a non-trivial topological structure, 
there being vacuum configurations of gauge fields (that is, pure gauge configurations) that 
cannot be continuously deformed into each other while keeping the system in vacuum. These 
configurations are classified into a discrete set, characterized by an integer 

n(A) = -|^y d 3 x e ^Tr(^A fe ), (2.12) 

where the temporal gauge Aq = is assumed. It is now straightforward to note that in 
transitions between topologically distinct vacua, B + L changes as A(B + L) = An. 

Under normal conditions, the different vacuum configurations are separated by a barrier. 
Processes taking the system from one vacuum to another are tunneling events and are de- 
scribed by instantons. However, these processes are highly improbable. At high temperatures 
the situation is different since thermal fluctuations can carry the system from vacuum to vac- 
uum. Moreover, the energy of the saddle point configuration of the barrier (called sphaleron) 
is proportional to the W mass and thus vanishes when the symmetry of the theory is restored. 
The baryon number violating processes are then unsuppressed. This can have severe conse- 
quences on the baryon number asymmetry in the universe and was one of the main motivations 
to study the electroweak phase diagram to high accuracy, as discussed in the introduction. 



2.3 Renormalization and notation 

In this thesis we will use dimensional regularization to regulate all the divergent integrals, i.e 
all the integrals are evaluated in d = 4 — 2e dimensions. Furthermore, we will employ MS 
renormalization scheme which amounts to writing the momentum integrals as 

f d d p 

J (2^F " M 

where A = /i (e 7 /47r)" 1 ^ 2 . All the couplings and observables all implicitly scaled to their 4d 
dimension by fi, so for example for pressure p(T) = ix 2e p{T) where [p(T) ] = 4 — 2e. 



A 2e [ -— 

47T 



d d p 
(2n) d 



(2.13) 
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The parameters of the theory will run with the renormalization scale. This is governed by 



, 2 (A) = , 2 (^) + ^(-^ 2 -^ 2 + 3 5y + 6A), 2 ln^ (2.14) 
A(A) = A(, ) + ^(^ + ^ + 3, V2 _^ A 

_^'2 A+12A 2_3 5 4 +6s 2^ ln A (215) 
2/a\ 2/\, 1 / 9 2 2 9 2 17 /2 \ 2 - ^ 



<£(A) = g^) + —^ Y -8gi--g'--g'^gi,ln-, (2.16) 

5 2 (A) = ^)-^<7 4 lA (2.17) 

5 ' 2 (A) = ff ,2 (/^) + ^ ,4 ln^, (2.18) 
<? 2 (A) = ff 2 (M) __I_ 5 4 ln A (219) 

Values of the parameters are fixed so that at /i = mz we have 

v 2 {m z ) = -m 2 ?, A(m z ) = -^G^mjj, 

g 2 Y {m z ) = 2V2G M m t 2 , <? 2 (m z ) = 4 V / 2G /1 m 2 ,, (2.20) 

</ 2 (m z ) = 4^2G M (m 2 - m 2 w ) , a s (m z ) = = 0.1187, 

where m# is the unknown mass of the Higgs boson, mw = 80.43 GeV, mz = 91.19 GeV 
and nit = 174.3 GeV are the masses of the W and Z bosons and the top quark, and 

= 1.664 • 10~ 5 GeV -2 is the Fermi coupling constant [52]. The matching of the val- 
ues of the parameters to the values of the physical observables is given here at tree level for 
simplicity. However, when precise numerical results are wanted (when determining the crit- 
ical temperature and the location of the endpoint of the first order phase transition line in 
chapter 5), the one loop relations given in [28] are used. We employ a power counting rule 
\ ~ g' 2 ~ g 2 ~ g%. ~ g 2 (unless otherwise stated) and assume the temperature to always be 
so high that the relation v 2 < g 2 T 2 applies. This power counting convention is then used to 
bring all the order of magnitude estimates to simple powers of g 2 . 

As will be discussed in the following chapter, we will use the imaginary time formalism 
to study finite temperature field theory. This in effect means that the spacetime will be 
Euclidean. We define the gamma matrices there so that 

{ 7/i ,7„} = 2<V, { 7 ^, 75 } = (2.21) 

and Tr 757 / i7^7p7< T oc e^ vpa + 0(e), where 5^, is the Kronecker delta symbol and t^ vpc is the 
totally antisymmetric tensor with €0123 = 1- 
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Chapter 3 

Field theories at finite temperature 



In this chapter we will briefly review the basics of finite temperature field theory. For a 
more detailed analysis, see for example [57]. We use here a notation such that tp refers to 
the collection of all the fields of a theory, but sometimes it is essential to make a distinction 
between bosons and fermions and we then use a notation such that <ft refers to all the bosons 
and ip to all the fermions. 

In the grand canonical ensemble, the equilibrium properties of any system at finite tem- 
perature are given by the density matrix, 



1 

— exp 



(3.1) 



where H is the Hamiltonian describing the dynamics of the system, TVj are all the conserved 
charges with fi-i being the corresponding chemical potentials and (5 = 1/T. In order to lighten 
the notation, we will from now on write Yl% ^i^i = V>N. Expectation values of physical 
variables are given as traces over the density matrix, (O) = Tr pO. 

In practice, however, it is more convenient to operate with the partition function Z(T, /Xj, V): 

Trp=l Z(T,m,V) = Tr exp [—(3 (H — fJ,N)] 

= J2^n\e^[-P(H-fiN)]\if n ). (3.2) 

Here ip n are the eigenstates of the operator H — fiN. Physical variables such as pressure p, 
energy E, entropy S and number of particles Ni are given by the standard thermodynamic re- 
lations. Defining thermodynamic potential (free energy) F by F(T, fii, V) = — T In Z(T, /Xj, V), 
we have 

dF dF 

P ~ ~dV 7 ~ ~df' 
3F 

Ni = -—, E = F + TS + J2^ N i- (3-3) 

i 

The partition function is most conveniently evaluated as a path integral. The matrix 
element (y n |exp[— (3(H — |txiV)]|<^ n ) can be regarded as a transition amplitude in imaginary 
time r = it, carrying states from r = to r = (3, when the "time evolution" is governed by 
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the operator H — jiN . Such a matrix element can be written as a path integral 



(Hexp [-/?(#- //TV)] bi> (3-4) 
= / VcpVir exp / dr / d 3 x (i7r(x, r)</?(x, r) — K(7r, <£>) + fjJ\f(ir, <p)) 

J ip(x.,0)=tpi UO J 

where 7r(x, r) are the canonically conjugate fields of </?(x, r) and H and are the Hamiltonian 
and particle number densities, respectively. Since in all the cases of interest to us, TC — fiM 
is at most quadratic in n, the integration over the momentum fields is Gaussian and can be 
performed immediately. The result is 

Z(T, n,V) = j Vip exp (J* dr J d 3 x £'(<p, <p)\ , (3.5) 

where the integration is over all the fields satisfying (anti-)periodic boundary conditions in r 
as discussed below. The Lagrangian jC'(ip, <fi) above is calculated from H — fjj\f and thus differs 
from the Lagrangian defining the vacuum theory when the chemical potentials are non-zero. 
More precisely, given H(ir, ip) — ^N(tt, ip), the associated Lagrangian will be 

C'(<p, if) = iir((p, ip)ip - H {ir(<p, if), if) + fjN (ir((p, <p), ip) , (3.6) 

where n((p,(p) is solved from 

ILP = ft? • (3 ' 7) 

Due to the trace in the definition of the partition function, the fields in the path integral 
must satisfy appropriate boundary conditions in the compact imaginary time. It follows that 
bosonic fields must be periodic, 0(x, 0) = 0(x, /?), while fermionic fields, being Grassmann 
variables, must be anti-periodic, Y>(x, 0) = — Y>(x, Hence it is convenient to expand the 
fields in Fourier series in imaginary time. From the boundary conditions above it follows that 

oo 

9 9(x,r)=T Y, ^"VnW, (3.8) 

n=— oo 

where the so called Matsubara frequencies are given by 

{Itl-kT for bosons, 
(3.9) 
(2n + 1)-kT for fermions. 

With this Fourier expansion it then becomes possible to reinterpret the finite temperature 
theory as a three dimensional field theory with infinitely many fields <£> n (x), n G N. 

The gauge fields, containing non-physical degrees of freedom, must be treated with care 
(for a review, see [17]). The gauge must be fixed in order to ensure that the integration is only 
over the configurations that are not gauge equivalent with respect to each other. Temporal 
gauge Aq = is a typical choice (though it is not a sufficient condition). Also, since the 
Gauss law is not included in the Hamiltonian equations of motion, it should be considered 
as a constraint on irf (the canonically conjugate fields of Af). In the absence of other fields, 
the Gauss law states that Dirf = 0. This constraint can be taken into account by means 
of Lagrange multiplier fields, for which we can impose periodic boundary conditions. In 
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the end we can identify the Lagrange multipliers as the temporal components of the gauge 
fields and after integrating over the momentum fields, the resulting Lagrangian will have 
the conventional covariant form of a gauge field theory. The remaining gauge freedom, due 
to invariance under gauge transformations that are periodic in the imaginary time interval 
r G [0, /?], can be removed by the Faddeev-Popov procedure, which introduces ghost fields 
to the theory. Although being Grassmann variables, they obey periodic boundary conditions 
since they are related to gauge transformations that are similarly periodic. 

As a concrete example, consider an SU(2) non-Abelian gauge field theory with a fermion 
and a scalar field transforming in the fundamental representation of the gauge group. The 
Lagrangian is given by (in Minkowski spacetime) 



(3.10) 



and the theory contains two linearly independent, mutually commuting conserved charges, the 
baryon number and the third component of isospin, 



B 

Qs 



d 3 x ip~/ i/j, 



d 3 x 



1 



(A^) 1 " r 3 $ - $V 3 D $ ) - € ibc A v > b G 



T fJ,U 



(3.11) 



Introducing chemical potentials [xb and hq for these charges and integrating over the canon- 
ically conjugate momentum fields we get for the partition function 



Z = I "Dtp exp 



J dr J d 3 x (C E + HBip7o4>) 



(3.12) 



where Ce is the Euclidean Lagrangian obtained from the original Lagrangian in Eq. (|3.10|) by 
going to imaginary time and by making a replacement [41,42,58] 



4 3 



4 3 



9 



(3.13) 



Gauge fixing and ghost terms are suppressed in the expression. We thus see that the chemical 
potentials related to the conserved gauge charges can be interpreted as background fields for 
the temporal components of the gauge fields. 

Although the relation between the temporal components of the gauge fields and the chemi- 
cal potentials related to the corresponding gauge charges can be obtained by a straightforward 
calculation, there is a simple argument that suggests it as well. Heuristically, consider a sys- 
tem which contains an external charge density q = Q/V = eN/V where e is the elementary 
charge and N is the number of external charged particles. This can be taken into account in 
the theory by adding a corresponding source term to the action: 



Z(T,q) 



Vip exp 



J dr J d 3 x (C E + iqA ) 



(3.14) 



The related chemical potential \xq will then be the conjugate variable to N, given by: 
T d 



Z(T,q) ON 
ie 



Z(T, 



(3.15) 

J drj d 3 x (£ E + iqA ) = -ie(A ), 
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suggesting that these chemical potentials can be associated with the expectation values of the 
temporal components of the gauge fields [58] . 

A condition for thermal equilibrium is that the free energy is stationary with respect to 
fluctuations around the expectation value of the gauge field (Aq). Thus, thermal equilibrium 
can be achieved only when the system is neutral with respect to gauge charges: 

dF dF 

= WT\~ir-~ Q - (3 - 16) 

d(A ) d(i Q 

While the chemical potential for the baryon number (in general for any charge whose conserva- 
tion is guaranteed by a global symmetry) can be chosen freely, we see that this is not the case 
for the chemical potentials for the gauge charges. Their values are determined by requiring 
the system to be neutral with respect to them. In general, their values will then depend on 
the chemical potential for the baryon number. Note that in the example above, the baryon 
number does not carry any gauge charge (i.e. for any the number of baryons with third 
component of isospin being +1/2 is equal to those with —1/2) and thus the system is neutral 
when fiQ = 0. This is not the case in general. In the electroweak theory, global charges carry 
a non-zero (hyper)electric charge and the chemical potential for the (hyper)electric charge will 
then be non-zero. 



3.1 Dimensional reduction 

Perturbative evaluation of the path integral in Eq. (j3.5|) at high temperatures is unreliable due 
to infrared divergences which arise when integrating over the bosonic zero modes </>o(x) [16,17]. 
In loop expansion these modes behave as massless particles in three dimensions which makes 
the infrared behavior of the integrals worse. On the other hand, integration over the non-static 
(i.e. u> n 7^ 0) modes is infrared safe, since temperature, providing these modes a mass, acts 
as an effective infrared regulator. Since there is a clear scale hierarchy between the static and 
non-static modes, it is natural to consider developing an effective field theory describing the 
static modes by perturbatively integrating out all the other modes. This is called dimensional 
reduction [23,24]. 

The dimensional reduction has many advantages. It essentially divides the problem of 
computing the path integral into two stages: 1. calculating the contribution from the non- 
static modes to the variable in question and constructing the effective theory, both tasks that 
can be carried out perturbatively, and 2. solving the effective theory, which one usually has to 
do with numerical computations. However, typically the effective theory is much better suited 
for numerical computations than the original four dimensional theory since it does not contain 
any fermions which are problematic to implement on a lattice. Also, since all the heavy scales 
have already been integrated out perturbatively, the lattice spacing can be taken to be larger 
(in physical units) which leads to larger lattice sizes and thus to a more precise treatment of 
the soft scales. 

The effective theory is constructed to be the most general theory describing the relevant 
degrees of freedom which has the required symmetries. In general, any non-renormalizable 
coupling will be suppressed by powers of the large scale that was integrated out, in this case 
7rT and thus it normally is sufficient to consider just the renormalizable effective theories. The 
parameters of the effective theory are mapped to those of the full theory by requiring that the 
effective theory reproduces all the static, bosonic Green's functions to the desired accuracy. 

For a generic gauge field theory, the effective theory resulting from integrating out the 
non-static modes will be a three dimensional gauge field theory with fundamental (Higgs) and 
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Figure 3.1: The ratio of the two mass scales in the electrostatic effective theory for an SU(2) + 
fundamental Higgs theory and for the standard model. The precise expressions for the masses 
are given in chapter and in appendix 1X1 



adjoint (temporal components of the gauge fields) scalars. The adjoint scalars correspond 
to the electrostatic modes of the original gauge field theory, E{ ~ F^i ~ diAo, and they 
obtain a thermal mass m 2 D ~ g 2 T 2 in this effective theory due to Debye screening of electric 
fields. Consequently, this theory is often called the electrostatic effective theory. Also the 
fundamental scalar mass acquires temperature dependent corrections. In theories where 
the fundamental scalar is responsible for a spontaneous symmetry breaking, such as the Higgs 
scalar in the standard model, these corrections are the cause of the symmetry restoration at 
high temperatures, m| ~ —v 2 + g 2 T 2 > when g 2 T 2 > v 2 . The three dimensional vector 
gauge boson, Ai, corresponding to the magnetic sector of the original theory (Bi ~ (-ijkFjk ~ 
€ijkdiAk) remains, however, massless, protected by gauge invariance. There nevertheless is a 
mass scale associated with it since in three dimensions the coupling becomes dimensionful, 

gl ~ g 2 T. 

At very high temperatures the adjoint and fundamental scalars are parametrically equally 
heavy, 7713 ~ ~ gT, and thus both of them should be considered on an equal footing. 
However, if the fundamental scalar is responsible for a spontaneous symmetry breaking, then as 
the temperature is lowered, it becomes increasingly light compared with the adjoint scalar, as 
shown in Fig. 13.11 This indicates the onset of the phase transition in which the symmetry gets 
spontaneously broken. Near the phase transition point we then have again a mass hierarchy 
with the adjoint scalars being heavy and it is natural to integrate them out, leaving an effective 
theory containing a fundamental scalar and three dimensional gauge fields. This theory can 
then be used to compute the characteristics of the phase transition. 

Even though the infrared problems related to screening of the electrostatic modes are tamed 
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by the Debye mass induced for the electric modes, also the electrostatic effective theory suffers 
from infrared problems in perturbation theory. These are related to screening of the magnetic 
fields. This can be studied by integrating out the massive adjoint and fundamental scalars 
to obtain an effective theory describing the magnetostatic degrees of freedom, namely the 
three dimensional gauge bosons. The resulting theory is confining and it cannot be solved 
with perturbative calculations. This places a limit on how far the perturbative expansion of 
different variables can be extended. 

We will consider the effect of introducing finite chemical potentials on the dimensional 
reduction in more detail in section 15.21 but let us briefly study some generic features of it. 
Chemical potentials for fermion numbers will change the form of the fermion propagators and 
thus the matching of fields and parameters of the different theories changes. This, however, will 
not change the structure of the effective theory. On the other hand, some of the symmetries 
of the four dimensional theory are explicitly broken by the finite chemical potentials, namely 
charge conjugation C, and this leads to new terms in the effective theory. More precisely, finite 
chemical potentials break C but preserve parity P and time reversal symmetry T and thus the 
theory is, in addition to being C breaking, also CP and CPT breaking (note that there is no 
controversy here since finite chemical potentials break Lorentz invariance). Depending on the 
field content of the system, the effective theory will then have a number of new terms. 
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Chapter 4 

Pressure of the standard model 



Pressure of a high temperature field theory is an important quantity both theoretically and 
phenomenologically. Defined by the path integral 

p(T) = lim ^ln [ VAVipVipV® exp(-S), (4.1) 

it is the most basic quantity that one can compute in a high temperature quantum field 
theory and is thus of special interest. From a phenomenological point of view, knowing the 
pressure (and the other thermodynamic functions that can be evaluated when we know the 
pressure) is important, for example, in cosmology. Expansion of the universe, governed by the 
Einstein equations, depends on the pressure and energy density of the matter it is composed 
of. Assuming the universe to be flat (which it according to observations is) and having no 
cosmological constant (contribution of the present day cosmological constant on the expansion 
is negligible in the early universe), the evolution of the scale factor R(t) is governed by 

^ 2 = f = "feiT), (4.2) 

| = -i^( e( T) + 3p(T)), (4.3) 

where e(T) is the energy density. Thus many theoretical predictions in cosmology depend on 
how well we know the pressure, energy and also entropy densities in the early universe. 

As an example, consider cold dark matter. Observations performed with the Wilkinson 
Microwave Anisotropy Probe (WMAP) suggest that as much as 23% of the matter in the 
universe is cold dark matter [59]. One possible candidate to explain cold dark matter is 
weakly interacting massive particles, or WIMPs, that might have been produced in the very 
early universe (for a recent review, see [60]). In order to be able to predict the density of 
WIMPs in the universe today, we need to know the equation of state, i.e. pressure, during 
the time when they decoupled from the thermal evolution of the universe [61]. It has been 
argued that differences of the order of 10% in the equation of state can alter the relic density of 
WIMPs as much as 1% [62]. It is possible that with the Planck satellite, designed to measure 
the power spectrum of the microwave background radiation with even higher precision than 
WMAP, we can measure the density of cold dark matter to this accuracy [63] and hence it is 
important to achieve the corresponding precision in theoretical calculations as well. 

In this chapter we will review the calculation of the pressure in the electroweak sector 
of the standard model, assuming temperature to be so high that the system resides in the 
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symmetric phase of the theory. The perturbative expansion is performed to order g 5 in the 
coupling constants, equivalent to three-loop order in loop expansion. In the next section, we 
analyze the structure of the expansion using a simpler SU(2) + Higgs theory, after which we 
consider the numerical results for both the simpler theory as well as for the full electroweak 
theory. 



4.1 Structure of the perturbative expansion of pressure 

For our purposes, full description of the electroweak theory requires five parameters: the gauge 
couplings g and g', the Higgs mass scale v 2 and self-coupling A and the top Yukawa coupling 
gy . In order to incorporate the QCD effects we have to add the strong coupling constant 
g s . The perturbative expansion of pressure becomes then a very complex function of all these 
parameters and to study the structure of the expansion is strenuous. Also, since the Higgs 
sector of the theory carries only 4 of the total of 106.75 degrees of freedom, the contribution 
of the Higgs to the pressure is not easily visible. 

In order to circumvent these problems, we will consider here a simpler theory that however 
has all the necessary properties, namely an SU(2) + fundamental Higgs theory defined by the 
Euclidean Lagrangian 

C = ^G^G^ + D^D^-vH^ + xf^^Y , (4.4) 

where G® v = d^A® - d v A^ + ge abc A b ^A c u and = - igr a A a ^/2. Characterized by 

only three couplings, the expansion for the pressure is easier to approach and the effects of 
the scalar sector are more evident. The expansion for the full electroweak theory is written 
down in appendix El an d the related numerical results are reviewed in section H~2l 

We will employ the framework of dimensional reduction, described in section 13.11 The 
pressure can then be written as (we implicitly assume that the limit V — > oo is taken) 

p(T) = pEW + ^lnJvAiDAoDQexpi-SE), (4.5) 

where Pe(T) is the strict perturbative expansion of the pressure of the 4d theory and Se 
defines the effective theory describing the n = bosonic Matsubara modes, 

= J d 3 x (jG%G% + \ {D.Alf + \m 2 D A a A a + \\ A (A a A a ) 2 

+Di&Di<5> + mf + A 3 (V$) 2 + h 3 A%A%&<^ , (4.6) 

where for the adjoint scalars Aq we have D^Aq = diA^ + g^e abc A\AQ. The matching of the 
parameters of this effective theory to those of the 4d theory has been done in detail to the 
desired accuracy in [2,28] and we will just review the results here. 

The contribution of the non-zero Matsubara modes (that is, the scale 2-kT) to the pressure 
is contained in the parameter pe(T) and in the couplings and masses of the effective theory 
Se- The goal is to evaluate the expansion of the pressure to order g 5 which determines the 
accuracy with which we need to know these parameters. For the couplings of the effective 
theory it is sufficient to use the tree level matching, 

g 2 = g 2 T, A 3 = AT, 

1 (4-7) 
h 3 = -g 2 T, \ A = 0{g% 
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Corrections of the order of g 4 to these would contribute to order g 6 in the pressure and thus 
we can neglect them now. The masses are needed to order g 4 : 
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Here we have explicitly separated the contribution coming from the Higgs scalar to the Debye 
mass mo by keeping the number of fundamental scalars ns general. This allows us to keep 
track of the effects of the Higgs sector. In the end one should set ns = 1. Order 0{e) terms 
are required for the masses since there will be terms of the form m 2 /e in the pressure of 
the effective theory. Note also that the mass of the adjoint scalar is renormalization group 
invariant to the order needed while the mass of the fundamental scalar develops a pole at 
order g\ in the effective theory that must be renormalized and therefore m\ runs. The mass 
counterterm, which can be obtained by the matching procedure, is 
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(4.10) 



The pressure pe(T) is obtained from the 4d theory by calculating its pressure in a strict 
perturbative expansion, i.e. evaluating the diagrams depicted in Fig. 14.11 This will lead to 
infrared divergences at order g A (at three loops) which are not canceled even after renormal- 
ization. They will be canceled only when the pressure of the effective theory is taken into 
account. Since we assume the temperature to be so high that v 2 < g 2 T 2 , we can expand the 
Higgs propagator in powers of v 2 . Such expansion is possible since Pe(T) counts contributions 
only from the non-static modes and integration over them is infrared safe. Thus the expansion 
of pe(T) in terms of v 2 is analytic. 
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Computing the diagrams in Fig. 14.11 (done in [2]), we get for Pe{T) 
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Here we have again explicitly written down the contribution from the Higgs scalar. The 
pressure p(T) is normalized so that p(0) = and the subtraction of the vacuum pressure, 
Pvac. ~ is taken into account in the ~ is 4 term of the expansion of p^,(T). Note that the 
expression in Eq. (|4.11j) for pe{T) is a high temperature expansion and thus is not suitable 
for studying the limit T — ► as such. The infrared divergences are manifest, corresponding 
to the terms ~ 1/e. 

Contribution from the soft scale gT to the pressure is contained in the pressure of the 
effective theory Se- Assuming that the adjoint and fundamental scalars are (parametrically) 
equally heavy (valid at very high temperatures), this is conveniently obtained by writing 

^ln J VA&AqV® exp(-S E ) = pm(T) + ^ In J VM exp(-5 M ) (4.12) 

where pm(T) is the contribution from the soft scales and the theory Sm describes the dynamics 
of the magnetic sector of the system, 
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(4.13) 



The theory is confining and cannot be solved using perturbative methods; however, it con- 
tributes to the pressure starting at order g and can therefore be neglected now. 

The pressure pm(T) is obtained by computing the loop expansion of Se to three loops. 
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Figure 4.1: Diagrams contributing to Pe(T) in the SU(2) + fundamental Higgs theory. The 
dashed lines correspond to the fundamental scalar and the wavy lines to gauge bosons. Dia- 
grams containing ghosts and counterterms are not drawn. 



The required diagrams are listed in Fig. 14.21 and the expansion is given by 
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The couplings of the effective theory, h 3 and A3, are substituted with those of the ori£ 
theory by using the matching given in Eq. (|4,7|) in order to make the expression shorter. 
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Figure 4.2: Diagrams contributing to pm(T) in the SU(2) + fundamental Higgs theory. The 
dashed lines correspond to the fundamental scalar, the wavy lines to the gauge bosons and 
the solid lines to the adjoint scalars. 



The pressure of the theory can now be written as the sum of the individual parts, 

p(T) = PE (T)+pu(T) + 0(g 6 ). (4.15) 

with pe(T) and pm(T) given by Eqs. (|4.11|) and (|4.14j) . respectively. Substituting the expres- 
sions for the masses mjj and ms to pm(T), we see that the poles in pe(T) and pm(T) cancel 
each other. However, a physical effect remains in the form of terms ~ g A \n(m/T) ~ g^lng 
where m refers to masses in the effective theory, 771,3 and rri£,. Additionally, the fact that the 
fundamental scalar mass runs in the effective theory leads to terms of the order g 5 In (A/7713) 
in the expansion. Dependence on the scale A vanishes order by order, as it should for a 
physical quantity, when running of the parameters is taken into account, but will not vanish 
completely unless an all-orders result is considered. The scale must therefore be fixed and we 
choose A = 2ttT. Sensitivity to changing the scale from this is studied later. 

4.1.1 Approaching the crossover transition 

Above we implicitly assumed that the fundamental and adjoint scalar masses are of the same 
order of magnitude (parametrically) . However, as the temperature of the system is decreased, 
the fundamental scalar becomes increasingly light compared to the adjoint scalar (see Fig. l3.1|) . 
Near the crossover transition in which the symmetry of the theory becomes spontaneously 
broken, the mass of the fundamental scalar will effectively vanish, m| ~ 0. This renders 
the computation of pm(T) unreliable near the crossover, as especially evidenced by the terms 
~ mjj/ms in the expansion. Thus pm{T) must be recalculated near the crossover. This is 
conveniently done by constructing a new effective theory. Since we have a scale hierarchy, 
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m 3 < m D, we can integrate out the adjoint scalars to obtain an effective theory for the 
fundamental scalars and 3d gauge fields, 
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(4.16) 



The couplings of this theory will to the required order be the same as the couplings of Se, 
93 = s|> -^3 = -^3- The fundamental scalar mass, on the other hand, does receive a correction, 
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We use a powercounting rule m 2 ~ fh 2 ~ g 3 T 2 (valid only near the crossover) and then 
the order g 4 corrections to the mass will contribute to the pressure at order g 5 ' 5 and will be 
neglected now. Note that fh 2 is renormalization group invariant to this order. 

The contribution from the effective theories to the pressure is now reorganized so that, 
instead of Eq. ()4.12j) . we have 



T 
V 
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= Pmi(T) + V wi(T) + ^ In J VAi exp(-S M ')- (4.18) 

The theory Sm' w iH differ from the theory Sm only by the matching of the parameters, the 
structure of the theories is the same. Its contribution to the pressure, of the order 0(<7 6 ), can 
again be neglected. 

The diagrams needed to calculate pmi(T) are given in Fig. 14.31 Those are computed in [3] 
and the result reads 
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where we have again substituted the matching of the couplings already to the expression. Note 
that, in addition to poles at order g 2 m 2 D that are canceled against the poles m.p^(T) and which 
correspond to infrared divergences related to screening of electric fields, there are also poles 
at order g 4 rri£) that are only canceled when the pressure PM2(T) is taken into account. These 
poles are related to the infrared divergences that appear in the theory given by Se when the 
mass of the fundamental scalar is taken to the limit 7713 — * 0. 

Since the fundamental scalar mass is small near the crossover, the pressure from Se> is 
only needed to two loops. The diagrams needed are shown in Fig. 14.41 and the result for the 
pressure reads 
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The poles here cancel against those coming from pe(T) and pmi(T). 

Collecting the results together, the pressure near the crossover is given by 



P(T) 



Pe(T) +PM1 (T)+p M2 (T)+T 4 0(g 
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(4.21) 
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Figure 4.3: Diagrams contributing to pmi(T). The solid lines correspond to the adjoint scalars, 
the dashed lines to the fundamental scalar and the wavy lines to the gauge bosons. 
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Figure 4.4: Diagrams contributing to PM2- The dashed lines correspond to the fundamental 
scalar and the wavy lines to the gauge bosons. 



This expression is well behaved when 777,3 — > as it should. Also, all the 1/e poles are canceled 
in the final result just like for the previous calculation and similar types of terms, of orders 
g 4 In g and g 5 In g are found in the expansion. Note that in this case, the presence of the terms 
~ g 5 In g is related to the infrared divergences that we encounter at the limit 7773 — > 0, not to 
the renormalization of the fundamental scalar mass in the effective theories. 

4.2 Numerical results 

In this section we will analyze the numerical consequences of the expansion derived in the 
previous section. We will start with the simpler SU(2) + Higgs theory, since the effects 
of the Higgs sector should be more evident in it. We will then turn to the full standard 
model for which the perturbative expansion is given in appendix [XJ The numerical values 
of all the relevant parameters are as given in section 12.31 with the Higgs mass chosen to be 
run = 130 GeV. The scale is chosen so that A = 2irT. 

The pressure of the SU(2) + Higgs theory is plotted in Fig. 14.51 at different orders of 
the perturbative expansion. Using two-loop effective potential calculations [18, 19] we have 
added the pressure of the broken phase to order g 3 to the picture to indicate where the phase 
transition (crossover) takes place. The result is normalized to the ideal gas pressure of massless 
particles po (T) , given by 



Po(T) 



7T 

~9 



(4.22) 



As can be seen from the figure, the pressure does not differ from the pressure of ideal gas 
by more than 2%. This reflects the fact that the theory is weakly coupled and thus the effect 
of interactions remains small. This is confirmed by the fact that the perturbative expansion 
seems to converge well: contribution of each new term in the expansion is smaller than that 
of the previous terms. 
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Figure 4.5: The pressure of the SU(2) + Higgs theory at different orders. 
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Figure 4.6: Scale dependence of the pressure of SU(2) + Higgs theory at different 
orders. Temperature is fixed to 500 GeV. 
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Figure 4.7: The equation of state parameter w(T) plotted for SU(2) + Higgs theory. The 
ideal gas result for radiation is given for reference. 



Convergence of the expansion can also be studied by considering the scale dependence 
of the result. Scale dependence should be reduced with each new order in the expansion if 
perturbation theory is valid. As can be seen in Fig. 14.61 this is partially so in this theory. 
Running the scale by 6 orders magnitude at a fixed temperature, the pressure changes only by 
about 2%. Also, scale dependence seems to be largest at order g 2 . However, the change in the 
pressure induced by running the scale is as big as is the general effect from the interactions. 
Also, the difference between the scale dependence at orders g 2 and g 5 is not large. In fact, the 
result is less dependent on the scale at orders g 3 and g 4 than at order g 5 . This is in accord 
with results obtained in [64] for the pressure of QCD. We can thus deduce that although the 
perturbative expansion seems to converge well, it can not be ruled out that higher order terms 
have as large a contribution to the pressure as the presently calculated terms. 

Since there is another explicit mass scale in the theory, v 2 , we can also expect the equation 
of state of this system to differ from that of massless particles, at least when the temperature is 
not too much above the scale v. This can be conveniently studied by computing the equation 
of state parameter w(T), defined by p(T) = w(T)e(T). For radiation w(T) = 1/3 while for 
non-relativistic matter we have w(T) = (pressureless dust) and for cosmological constant 
w(T) = — 1. In Fig. 14.71 we have plotted w(T) for matter described by the theory in question. 
As can be seen, for high temperatures the system behaves very much as radiation, but as 
the temperature decreases, we see deviation from w(T) = 1/3. This is to be expected since 
the terms v 2 T 2 and v A become increasingly important as temperature is lowered, making the 
pressure to deviate from the form p ~ T 4 . The difference to w(T) = 1/3 of radiation is, 
however, still small. Although w(T) seems to grow near the crossover transition, one expects 
it to decrease soon after the transition since the system will contain a number of massive 
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Figure 4.8: The pressure of SU(2) + Higgs theory as plotted all the way down to the crossover 
temperatures. Unphysical singularity in the high temperature result put(T) is manifest while 
the consistent calculation taking into account the lightness of the fundamental scalar, denoted 
by ppt(T) in the graph, is seen to behave well. 



fields whose contribution to the pressure is negligible compared with their contribution to the 
energy density, p(T)/e(T) = T/m, m 3> T. 

As we decrease the temperature, we eventually reach the crossover transition and the 
expression used for pressure at high temperatures (denoted now as pht(T)) ceases to be valid. 
This is seen in Fig. 14. 81 in which we see an unphysical singularity for the pressure when plotted 
using the high temperature expression given in Eq. (|4.15|) . The correct result for the pressure 
in that region, j>pt(T), defined by Eq. 1)4.21)1 . is given by the solid curve in Fig. 14.81 However, 
the singular behavior of pwv{T) is isolated to a narrow range of temperatures. This can be 
understood by considering the leading order terms responsible for the singularity, the terms 
of the form ~ rn^/m^ in the expansion of pressure. We get for the singular terms 

Psingular(T) _ 135 fl 6 

p (T) 4096vr5y3^7^V5T' 1 ' > 

where To is the temperature such that m^{To) = and 5T is the deviation from that, 5T = 
T — To. Since the numerical factor of this term is very small, the temperature must be very 
close to To for any effects from this term to be manifest. 

We can also see that ppt(T) runs consistently somewhat below the high temperature result. 
This is due to the fact that when inserting the mass mj from Eq. 1)4. 17|) to the expansion of 
the pressure PM2(T) in Eq. 1)4.20)) . we effectively resum the class of diagrams in Fig. 14.91 
to the pressure. The high temperature result, instead, would correspond to expanding the 
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Figure 4.9: The type of diagrams that are resummed to the pressure ppt- The dashed line 
corresponds to the fundamental scalar and the solid lines to the adjoint scalars. 



mass m| in powers of h^mn/m^, the leading order given by m| = m§, when inserted to the 
expansion of pressure. At high temperatures, this expansion can be carried out since then 
h^mi) /m| ~ g <C 1, but near the crossover hsmD/m^ ~ 1 according to our power counting 
rules and such expansion cannot be performed then. 

Let us next consider the pressure of the full standard model. The expansion in powers of 
the coupling constants is given in appendix^ The pressure, including the QCD contribution, 
is plotted at different orders of perturbation theory in Fig. 14.101 normalized again to the 
pressure of ideal gas of massless particles, 

p (T) = 106.75^T 4 . (4.24) 

It is seen to significantly deviate from the ideal gas result, by up to 15%. The convergence 
of the expansion is also much worse than that of the SU(2) + Higgs theory. This was to 
be expected since there are two large couplings in the theory, the strong coupling constant 
9s{ m z) ~ 1-5 and the top Yukawa coupling g\{mz) ~ 1-4. The scale dependence of the result 
is not plotted here since one does not expect it to differ from the QCD results [64,65]: QCD 
degrees of freedom constitute the majority of all the degrees of freedom in the standard model 
and, moreover, the scale dependence coming from running of g s has by far the largest effect. 

The pressure near the electroweak crossover is plotted in Fig. 14.111 similarly normalized. 
The unphysical singularity in the pressure when plotted using the expression valid at high 
temperatures is again manifest, but the calculation that takes into account that the funda- 
mental scalar is light around the crossover behaves smoothly. Note that the singular behavior 
is again isolated to a very narrow temperature range, even more so than in the SU(2) + Higgs 
theory. The reason is that this effect stems from the fundamental scalar sector which, as 
already noted, carries just a tiny fraction of all the degrees of freedom and thus the numerical 
factor in front of the singular terms (normalized to the ideal gas pressure) is even smaller (see 
Eq. 

We can study the effects that the scalar sector has on the pressure by altering the mass 
of the Higgs boson. Plotted in Fig. 14.121 we see that changing the Higgs mass from 130 GeV 
to 200 GeV has a minimal effect on the pressure, again reflecting the small relative number 
of degrees of freedom. It is then not necessary to repeat the numerical analysis here for a 
number of different Higgs masses. 

Using the result for the pressure, we can compute also other thermodynamic variables. 
In Fig. 14.131 we have plotted the effective number of bosonic degrees of freedom in terms of 
pressure (/ e ff)) energy density (g c g) and entropy density (h e g). Those are defined by 

7T 2 7T 2 27T 2 

P(T) = f cS ^T\ e(T) = g*T\ s(T) = h^T* (4.25) 
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Figure 4.10: Pressure of the standard model at different orders of perturbation theory. 
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Figure 4.11: Pressure of the standard model near the electroweak crossover. 



and for ideal gas their values would be equal, f*g = g 1 ^ = h 1 ^ = 106.75. As can be seen, the 
effect of interactions on the energy and entropy densities is comparable to that they have on 
the pressure. 

Finally, in Fig. 14.141 we plot the equation of state parameter w(T) and the speed of sound, 
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Figure 4.12: Relative difference between pressures at iuh = 130 Gev and mjj = 200 GeV. 
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Figure 4.13: The effective number of degrees of freedom in terms of pressure, energy density 
and entropy density, as denned in the text. For ideal gas each of these would be equal to 
106.75. 



(? S (T) = p'(T)/e'(T) (where ' refers to derivative with respect to temperature). They both 
equal 1/3 for ideal gas. Unlike for pressure, energy density and entropy density, the effect of 
interactions on these variables is observed to be rather weak. We do see deviation from the 
ideal gas results, but the effect is not much more than 2%. We can thus deduce that though 
the thermodynamic potentials do significantly deviate from the ideal gas results, the matter 
still seems to behave very much like radiation. 

In this thesis, the computation of the pressure is carried out assuming the system resides 
in the symmetric phase of the theory. However, it is of phenomenological interest to evaluate 
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the pressure also in the broken symmetry phase. This can be done by extending the effective 
potential calculations to the corresponding precision, but it is a highly non-trivial task due to 
the complicated structure of the theory in the broken phase. Phenomenological estimates for 
the pressure at temperatures between the electroweak scale and the QCD scale have, however, 
been presented [66]. 
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Chapter 5 



The electroweak phase diagram at 
finite chemical potentials 

As discussed in the introduction, much of the research on electroweak thermodynamics has 
concentrated on studying the properties of the electroweak phase transition. In this chapter, 
we will review the effect of finite chemical potentials related to conserved fermion numbers on 
the phase diagram. The motivation is to understand the phase structure of the electroweak 
theory in more general terms, but the computation may have phenomenological interest as well. 
Although the baryon number asymmetry in the universe can be estimated with observations 
(baryon-to-photon ratio is of the order of 10 -10 ), measuring the lepton number asymmetry 
in the universe is more difficult since neutrinos interact only weakly with other particles. 
Consequently, a large lepton number asymmetry, residing in neutrinos, has presently not 
been ruled out. For this reason, we will in the next section briefly review some aspects of 
neutrino cosmology. We then move on to consider the computation of the phase diagram of 
the electroweak theory in presence of finite chemical potentials. 

5.1 Lepton asymmetry in the universe 

Neutrinos interact with other particles only via weak interactions and for this reason they 
decouple from thermal evolution of the rest of the universe at an early stage. Assuming 
the neutrinos are massless, the interaction rate of processes keeping the neutrinos in thermal 
contact with other particles (such as ve <-> ve and vv <-> ee) is roughly given by T = na ~ G^T 5 
where n is the number density of neutrinos and we have estimated the cross section a to behave 
as a ~ G^T 2 . This rate is of the same order as the expansion rate of the universe when the 
temperature is about 1 MeV, assuming that there is no large asymmetry between the numbers 
of neutrinos and antineutrinos. After decoupling, neutrinos retain a thermal distribution, but 
the temperature T u of the neutrinos, decreasing as the universe expands according to T v ~ 1/R, 
may differ from the temperature of the rest of the universe (photons). Indeed, shortly after the 
neutrino decoupling, the temperature of the universe drops below the electron mass threshold 
and entropy carried by the electron-positron pairs is transferred to photons. Consequently, the 
temperature of photons will be boosted by a factor of (11/4) 1 / 3 compared to the temperature 
of the neutrino background. Since the temperature of the microwave background radiation 
today is about 2.735 K, we get that the corresponding temperature of the neutrino background 
is about 1.95 K. This corresponds to a neutrino background of about n v ps n p pt 56 neutrinos 
per cm 3 in the universe today. 
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However, there are no direct observations concerning the neutrino degeneracy, i.e. asym- 
metry between neutrinos and antineutrinos, in the universe. Upper limits for the neutrino 
degeneracy can be obtained by considering big bang nucleosynthesis (BBN) which is sensitive 
to a degeneracy in electron neutrinos v e . Excess of v e with respect to v e will induce changes 
in the beta reactions leading to lower neutron to proton ratio on which the abundances of 
primordially produced light elements depend. This constrains the chemical potential of v e to 
lie between —0.01 < n Ve /T Ve < 0.22, where T Ve is the temperature of electron neutrino back- 
ground. Asymmetry in the [i and r neutrinos affects the BBN by hastening the expansion of 
the universe and this bounds the corresponding chemical potentials by T /T v T | < 2.6 [67]. 
One can also obtain limits for the neutrino degeneracy by considering the power spectrum of 
the microwave background radiation. Based on the data from WMAP, Lattanzi, Ruffini and 
Vereshchagin [68] find that the neutrino degeneracy is constrained tobeO<|/x/T|<l.l and 
claim that statistical fits, in fact, prefer a large neutrino degeneracy, fi/T ~ 0.6 

There are a number of models that could explain a presence of a large lepton number 
asymmetry in the universe without leading to a comparable baryon number asymmetry, see 
for example [69,70]. An obvious problem is to avoid the transformation of any generated net 
lepton number to a comparable net baryon number via sphaleron mediated processes that are 
unsuppressed when the electroweak symmetry is unbroken. The common argument is that a 
large lepton number asymmetry would prevent the restoration of the electroweak symmetry 
even at high temperatures [36] and thus the sphaleron processes would in fact be suppressed. 
One can also consider a situation in which there are large net lepton numbers Lj within each 
family, but in such a way that the sum of the lepton numbers is small, L e + L^ + L T <C L e ^ T . 
This requires fine tuning and thus is not natural. 



5.2 Dimensional reduction at finite chemical potentials 

As discussed previously, the exactly conserved global charges in the standard model are the 
baryon - lepton number charges, 

Xi = - Li, i = e,fi,T, (5.1) 
where i refers to different families and B and Lj are defined by 

B = \Yj J d3x Qc,i 70 q c ,i, (5.2) 

c,i 

Li = j d 3 x (eryo ej + ^70 a L v%) , ax, = ^(1 - 75). (5.3) 

Here c refers to color, q c ^ are the quarks and and Vi are the electron type lepton and the 
corresponding neutrino of each family. We assign chemical potentials ^ to these conserved 
charges, but for computational reasons, it is convenient to introduce separate chemical poten- 
tials for the baryon number B and for the lepton numbers Lj, namely \ib and //j^. These are 
then related so that 

f*B = g VLi = -M- (5.4) 

i=e,fi,r 

For this reason, we will refer to the chemical potentials /ij as "leptonic chemical potentials" . 
Note that each colored quark carries a baryonic chemical potential /j,b/3- 
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Dimensionality Terms 

GeV 1/2 Wo 

GeV 3/2 i&A%T a <S>, i&B $, iB A%A% 

GeV 2 e ijk BiF jk , e ijk (a«G% - \ge ahc A a l A)Al 



Table 5.1: The dimensionally lowest order terms violating CP and CPT invariances in the 
electroweak theory. 



Additionally, there are conserved gauge charges, but as discussed in chapter |31 the cor- 
responding chemical potentials can be absorbed to the temporal components of the gauge 
fields and are thus not manifest in the path integral [41,42,58]. The equilibrium properties 
of the standard model in the presence of the chemical potentials \X{ are then described by the 
partition function 
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H B B + 



(5.5) 



where S is the Euclidean action of the electroweak theory. 

We are now interested in computing the phase diagram of the electroweak theory when 
the chemical potentials m are non-zero. A proper approach at high temperatures is again the 
framework provided by dimensional reduction. As discussed in the previous sections, close 
to the phase transition the relevant light degrees of freedom are the fundamental scalar <3? 
and the magnetostatic gauge bosons A?, Bi. Construction of the effective theory describing 
these fields is performed in detail in [1] and we will merely review the generic features of the 
computation here. 

As pointed out in section 13.11 introducing chemical potentials to the system will reduce 
the number of discrete symmetries the theory has (or at least may have) otherwise. Namely, 
the theory becomes CP and CPT breaking. The effective theory may then contain terms 
that break these symmetries but which still respect three dimensional gauge- and rotational 
invariances (for a review on how different fields transform under the discrete symmetries, 
see [71]). The most general renormalizable theory describing the three dimensional gauge 
fields and the fundamental scalar in question can then be written as 



+ae ijk (A1G% - ^g 3 e abc A* A)A%\ + /•),.. 



(5.6) 



Here the Chern-Simons terms (terms proportional to 5 and a') break CP and CPT and thus 
vanish automatically when chemical potentials are not present. However, they appear in the 
effective theory if we have a finite chemical potential for the non-conserved B + L [72] . Since 
we consider only the strictly conserved B — L numbers, the Chern-Simons terms will not be 
present in the effective theory, i.e. a = a' = 0. The structure of the effective theory is thus 
the same as it is when the chemical potentials vanish, only the matching of the parameters 
53> 9'ii ^3 an d A3 to the physical variables changes. The matching is given in [1]. 

When constructing the theory describing the light degrees of freedom, we again need to 
first develop an effective theory containing the electrostatic modes Aq and Bq as well. Unlike 
the theory describing the light degrees of freedom, the structure of the electrostatic effective 
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theory does change when chemical potentials are introduced. The dimensionally lowest order 
terms (that are not present when \i{ = 0) that must be taken into account in this theory are 
listed in table 15. ll for the standard model. The Lagrangian of the electrostatic effective theory 
is then 

£ E = + ~^ F ij F ij + (A*)*(A$) + ml&ff* + A 3 (<I> t <i>) 2 

+ \{D % Atf + \mlA a A a + \\ A {AlAtf + i(d tJ B ) 2 + \m%B B, 

+h 3 &$A a A a Q + ti 3 &$B B - ^g 3 g' 3 B &A%T a $ 
+ Kl B + k 3 B$ + p&A%T a $ + p'&$Bo + pgBqA^AI 
+ae ljk (A?G* k - ^e^A^ + a'e^B^ 



(5.7) 



The coefficients for the Chern-Simons terms vanish identically again. Of special importance 
is the term k\Bq. Having such linear terms in the Lagrangian leads, in equilibrium, to con- 
densates of the corresponding field, to leading order 



(Bo) 



4^0. 



m 



(5.8) 



As discussed before, this implies a finite chemical potential for the hypercharge. The linear 
term k±Bq in the Lagrangian thus ensures neutrality of the system with respect to hypercharge. 
Performing the matching (see [1] for details) we get to leading order for Kf. 



2f2 



TT 



ttT V 9 7T 2 T 2 



(5.9) 



Matching of all the other terms to physical parameters is given to high accuracy in [1]. 



5.3 The phase diagram 

The electroweak phase transition is described by the effective theory £ c g. (with 5 = a 1 = 0) 
introduced in the previous section. Due to infrared divergences encountered in perturbative 
calculations, this theory must be solved with numerical computations. To this end, lattice 
Monte Carlo studies have been carried out [29-33] and the phase diagram of the theory has 
been reliably solved. 1 In practice, one can neglect the U(l) subgroup and just consider SU(2) 
+ Higgs theory [31]. 

It is convenient to introduce dimensionless quantities 

% 9t 9s 

and let the coupling g| to provide the dimensions. Of these, z is essentially constant, z ~ 0.3 
(= for SU(2) + Higgs theory), and y determines the phase of the system: at tree level, if 
y > the system resides in the symmetric phase while if y < the symmetry is spontaneously 
broken. Thus, at tree level the critical temperature is determined by setting y = 0. The 
remaining variable, x, depending on the value of the Higgs mass, parametrizes the theory. 

1 We can apply the same numerical computations now, since the structure of the effective theory remains 
the same when finite chemical potentials are introduced. 
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Figure 5.1: The phase diagram of the effective theory in terms of the dimensionless variables 
x and y. The dashed line corresponds to the perturbative result based on two-loop effective 
potential calculations. 



The phase diagram of the theory in terms of the dimensionless variables x and y is plotted 
in Fig. 15.11 Both the perturbative result (based on two-loop effective potential calculations 
in the effective theory, [73]) and numerical results are presented. The perturbative result and 
the numerical results agree well for small values of x (which amounts to small Higgs masses 
and chemical potentials, as will be discusses shortly), but while perturbative computations 
suggest that there is a first order phase transition for all values of x, numerical studies show 
that there is an endpoint to the critical curve at x ~ 0.1 above which there is only a crossover 
transition. 

In order to map the phase diagram in Fig. l5.ll to a phase diagram in terms of temperature, 
chemical potentials and the Higgs mass, we need to apply the matching between the parameters 
of the effective theory and the physical 4d theory, given in detail in [1]. The critical temperature 
T c = T c (fj,i,m H ) is then given by y(T c , m,m H ) = y c (x(T c , m, tuh)) using y c (x) in Fig. EHJ In 
order to simplify the study, we assume all the leptonic chemical potentials to be equal to each 
other, fa = fi\/i, i = e, fi, r. Then, to leading order 



These expressions can be used to qualitatively understand the effect of the chemical potentials 
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Figure 5.2: The phase diagram of the electroweak theory in terms of the Higgs mass and 
temperature. 



on the phase diagram. As can be seen, the chemical potentials lead to decreasing the value 
of y and increasing the value of x. Thus, chemical potentials tend to break the symmetry of 
the theory and, consequently, the critical temperature will increase. This is in accord with 
results already obtained by Linde [36]. Moreover, the increasing value of x implies that the 
transition becomes weaker and that for sufficiently large chemical potentials there will be no 
phase transition at all, regardless of the value of the Higgs mass. 

These qualitative considerations are verified by the more accurate treatment of matching 
the phase diagram of the effective theory in terms of physical parameters. The electroweak 
phase diagram in terms of the Higgs mass and temperature is plotted in Fig. l5.2l for a number 
of different values of the chemical potential [i. The location of the endpoint of the first order 
phase transition line is seen to move to smaller values of the Higgs mass as the chemical 
potentials are increased, indicating that the chemical potentials make the transition weaker. 
At the same time, the critical temperature is slightly increased. Note also that, for some 
values of fj,, there is a first order phase transition only if the Higgs mass is within some 
range m^ wer limit < mu < m^ pcr limit . The reason is that, for a fixed and sufficiently large 
value of the chemical potential, the value of x along the first order phase transition line 2 , 

2 Of course, there is no first order phase transition line if x > Sendpoint ~ O.f. We can, however, formally 
continue the curve fitted to the lattice results in Fig. 15. f I to values x > Xendpoint and we refer to that continuation 
as the critical curve in that region of the parameter space. One should not attach any physical meaning to 
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Figure 5.3: The value of x c = x(T c ,^,m#) as a function of run for a number of different 
values of the chemical potential [i. 





endpoint 

m H 


rpendpoint 


GeV 


72 GeV 


109 GeV 


15 GeV 


71 GeV 


108 GeV 


30 GeV 


66 GeV 
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Table 5.2: The location of the endpoint of the first order phase transition line for a number 
of different values for the chemical potential. 



x c = x(T c , mn), is actually a decreasing function of mn for small values of ran , see Fig. 15.31 
Thus, although x(T c , /i, mn) > ^endpoint ~ 0.1 for small values of m# when /j, is large enough, 
it may decrease so much as mn is increased that it becomes smaller than x ent ipoint f° r some 
values of ma- 
in Fig. I5.4l we plot the phase diagram in terms of the chemical potential and temperature. 
As can be seen, there is again an endpoint to the first order phase transition lines, regardless 
of the value of the Higgs mass. If the chemical potential has an absolute value of fi > 50 GeV, 
there is no phase transition in the electroweak theory for any value of the Higgs mass, just 
a crossover transition. Finally, in table 15.21 we list the values of the Higgs mass and critical 
temperature corresponding to the endpoint of the first order phase transition line for a number 
of different values of the chemical potential. 

There is a simple physical interpretation of the results. Introducing a finite chemical 
potential for the conserved fermion numbers will in general lead to a non-zero (hyper)charge 

such a continuation, it merely serves as a convenient tool for analysis. 
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Figure 5.4: The phase diagram of the electroweak theory in terms of the chemical potential 
and temperature. 



density in the system. However, the system must be neutral with respect to gauge charges in 
equilibrium and thus a finite chemical potential for the hypercharge is generated to guarantee 
this. This effectively induces a finite chemical potential for the scalar doublet as well since 
it carries a hypercharge. As well known, a finite chemical potential for bosons leads to Bose- 
Einstein condensation 3 and thus, in this case, to symmetry breaking. To overcome this, even 
larger temperatures are needed to restore the symmetry of the theory. 

The computation of the phase diagram is by construction limited to small chemical po- 
tentials. The reason is that, in principle, dimensional reduction requires there to be a large 
scale hierarchy and the temperature to be the largest scale in the system. If there are other 
comparable scales, then the effective theories may contain non-renormalizable terms that are 
not suppressed by any large scale hierarchy and which therefore cannot be excluded from the 
theory. This would spoil dimensional reduction. 

Within QCD, it has been studied in more detail how large chemical potentials one can 
consider in this framework. Perhaps surprisingly, it turns out that one can allow the chemical 
potentials to be as large as \i < 4T [74]. Vuorinen [75], computing the QCD pressure and 
comparing calculations in the limit /i/T < 1 and at T = 0, also concludes that dimensional 
reduction seems to work well for very large values of n/T. Parametrically, dimensional reduc- 

3 For relativistic bosons, condensation happens when = m, where /x is the chemical potential and m is 
the mass of the bosons. For non-relativistic bosons \x < and condensation happens in the limit fi — > 0. 
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m H = 20 GeV 
m H = 38GeV 
m H = 56 GeV 
m H = 68 GeV 



tion works for T <C \x if the mass scales of the effective theories are smaller than the scales 
that were integrated out, ir 2 T 2 » m 2 D ~ g 2 s {T 2 + fi 2 /ir 2 ) ~ g 2 y} /it 2 , i.e., for fi/T <C TT 2 /g s - 
However, whether or not similar conclusions can be drawn within electroweak physics is an 
open question. 
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Chapter 6 



Comparison of electroweak and 
QCD thermodynamics 



In this thesis we have concentrated on electroweak thermodynamics. However, also finite 
temperature QCD is relevant, both from theoretical and phenomenological viewpoints: with 
ongoing experiments at RHIC in Brookhaven (and in near future at LHC at CERN) we are in 
a position to create and study strongly interacting matter [76]. In this chapter, we will briefly 
consider the similarities and differences between the QCD and electroweak thermodynamics. 

6.1 Comparison of the pressures 

Perturbative evaluation of the pressure of QCD with massless quarks has a long history [64, 
65, 77-83] and a number of theoretical breakthroughs has been achieved while extending the 
expansion further (for a brief review on the history of this topic, see [84]). Today, the expansion 
is known to the last perturbatively calculable term [75,82]. The generic structure of the 
expansion, as far as perturbative computations can extend, is 



Computation of the order g^ term cannot anymore be carried out within perturbation theory 
[16, 17] and to obtain it requires numerical input. First steps in computing it have already 
been taken [85]. 

Some general observations can now be made of the expansion. Although, in principle, 
there is a mass scale inherent to QCD, the Landau pole Aqcd> 



it only arises in quantum theory. At tree level the only mass scale is provided by temperature. 
Hence, the expression for the perturbative expansion of pressure deviates from the form p(T) ~ 
T 4 only by logarithmic terms. This should be compared to the expression for the pressure 
in the full standard model. There we have another explicit mass scale at tree level and, 
consequently, the expansion of the pressure has a different polynomial structure altogether, 
to leading order p(T) ~ T 4 + v 2 T 2 + v A + 0(v®). The complete expansion has an even richer 
structure, containing roots and logarithms of linear combinations of v 2 and T 2 , as seen in 
chapter ^ Similar structure applies to the couplings as well: in QCD with massless quarks 



PQCd(T) 



co + c 2 g 2 s + (c 4 + c' A In 5 S ) g 4 + c b gl + c 6 g® lng s + O (gf) . 



(6.1) 




(to one-loop) 



(6.2) 
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there is only one coupling constant, g s , and thus the expansion is a polynomial in powers and 
logarithms of g s . In the electroweak theory there are multiple couplings, leading to non-trivial 
combinations of them and thus the expansion is not a simple polynomial in the couplings and 
their logarithms. As a more subtle difference, there are no terms of the order In g s in the 
expansion of QCD pressure. This is related to the fact that the emerging effective theories are 
finite to the order needed and thus require no renormalization. In the electroweak case, the 
fundamental scalar mass runs in the effective theories at order g^ and this leads to terms of 
the order g 5 In g in the expansion. 

The perturbative expansion of QCD pressure converges poorly due to the magnitude of 
the strong coupling. Especially, evaluating the pressure near the QCD phase transition is 
in principle impossible by using perturbative methods since the critical temperature of the 
transition is close to the Landau pole of the coupling. The expansion of the electroweak 
pressure behaves much better in this sense: the theory is weakly coupled and the expansion 
appears to converge well. Moreover, electroweak pressure can be studied with perturbative 
computations even near the electroweak phase transition since the transition (crossover) is 
driven by the Higgs scalar and not by the theory becoming confining. However, as was seen 
previously, even in the electroweak case it cannot be ruled out that higher order corrections 
have a comparable contribution to the pressure. 



6.2 The phase diagrams 

Both QCD and the electroweak theory contain a phase transition. The exact properties and 
critical temperatures of the transitions depend on the chemical potentials and the values of 
parameters of the theories. In this section, we will compare the phase diagrams of these 
theories. The electroweak phase diagram is considered in terms of the leptonic chemical 
potentials and the theory is parametrized by the Higgs mass, while the QCD phase diagram 
[86-88] is considered in terms of the baryonic chemical potential and the theory is parametrized 
in terms of the strange quark mass (i.e. whether the system behaves as Nf = 2 of as Nf = 3). 

Consider first the cases when the masses parameterizing the theories are small. Then in 
both theories there is a first order phase transition, see Fig. 16.11 However, as the chemical 
potentials are increased, the critical temperature of the electroweak phase transition increases, 
while the critical temperature of the QCD phase transition decreases. Thus the responses of 
the systems on introducing the chemical potentials are opposite. 

Some interesting thermodynamics can be deduced from this. The Clausius-Clapeyron 
relation states that 

dT m - n 2 , q\ 

~T = > (6-3) 

where all the quantities are measured along the phase transition line and rii and Si are the par- 
ticle number and entropy densities in each phase, respectively. Since the entropy is maximized 
in the high temperature phase, we can deduce from the relation above that in an electroweak 
system, since dT/dfi > 0, the lepton number density in the low temperature phase (the broken 
symmetry phase) is higher than in the high temperature phase. In QCD the opposite is true, 
baryon density of the high temperature phase is higher. 

Another feature of the phase transitions that we can see from Fig. 16.11 is the response to 
increasing the masses that parametrize the theories. This is also depicted in Fig. 16.21 There 
is no first order electroweak phase transition if the Higgs mass is big enough, regardless of the 
value of the leptonic chemical potential. Also, if the leptonic chemical potential is larger than 
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Symmetric phase 



m H -m w 




m H « m w 



Broken symmetry phase 



Quark-gluon plasma phase 



m » m 




Hadronic 



Figure 6.1: Schematic plots of the electroweak (on the left) and QCD (right) phase diagrams 
in terms of temperature and the relevant chemical potentials. The solid lines correspond to 
the critical lines for a number of different Higgs/strange quark masses and the dotted lines 
indicate the location of the endpoint of the first order phase transition line as the masses are 
varied. The arrows indicate the order in which masses increase along the dotted lines. The 
low-T/high-^ phases of QCD (and also of electroweak theory), such as color superconducting 
phases, are not included in the picture. 



some critical value (fJ-L,c ^ 50 GeV), > f^L,c, there is no first order transition for any value 
of the Higgs mass. On the other hand, in QCD, if the baryonic chemical potential is large 
enough, fis > Hb,c, there is always a first order phase transition even in the limit m s — > oo. 
Note, though, that if \ib is very high and T is exponentially small, T ~ g~ 5 exp(— l/g s ), the 
low temperature phase of the system is not the normal hadronic phase but something more 
exotic (color superconductivity). Similar phenomena may be present in an electroweak system 
as well. In the present considerations, we assume fiB (a*l) to be small enough so that the low 
temperature phase is the hadronic phase (normal broken symmetry phase). 
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Figure 6.2: Schematic plots of the electroweak (on the left) and QCD (right) phase diagrams 
in terms of temperature and the masses parameterizing the theories for a number of chemical 
potentials and hb- The arrows indicate the order in which the chemical potentials increase. 



46 



Chapter 7 

Conclusions 



In this thesis we have studied some aspects of electroweak thermodynamics, namely the pres- 
sure of electroweak matter and the structure of the electroweak phase diagram. Due to lack of 
any experimental input about the properties of such matter, the motivation has been to un- 
derstand the theory better. The electroweak theory, as a part of the standard model, describes 
nature to high accuracy and thus it is important that we understand the theory completely. 

In this thesis, it was found that the pressure of matter composed of elementary particles at 
very high temperatures (T > 100 GeV) deviates significantly from the ideal gas pressure that 
is commonly used, for example, in cosmological computations. However, a great part of the 
deviation is due to the strong interactions and the contribution from the purely electroweak 
degrees of freedom is seen to be close to the ideal gas estimate. From a theoretical point of 
view, the perturbative expansion of the electroweak pressure exhibits, however, some novel 
features that were not encountered in previous computations of pressure of gauge field theories. 

Concerning the phase diagram of the electroweak theory, it was found that introducing 
finite chemical potentials for the conserved particle numbers leads to a weaker transition (when 
the parameters of the theory, i.e., Higgs mass, are such that a first order phase transition 
is possible) and that for large enough chemical potentials there is no phase transition at 
all, regardless of the value of the Higgs mass, just a crossover. To achieve this result, a 
generalization of dimensional reduction methods to finite chemical potentials was needed in the 
context of electroweak interactions. The resulting theories can be used to study the properties 
of electroweak matter at high temperatures and finite chemical potentials, in general. 

There is still much to learn about electroweak thermodynamics. The pressure of the stan- 
dard model in the broken symmetry phase with all the particles taken into account with their 
proper masses has not been evaluated to high precision, although it is of phenomenological 
interest in cosmology. The structure of the electroweak phase diagram at large chemical poten- 
tials and (exponentially) low temperatures is fairly unknown as well. Such considerations have 
been carried out within QCD and they suggest that the phase structure might be very exotic. 
It is therefore to be expected that research within electroweak thermodynamics remains active 
in the future as well. 

The author would like to thank M. Vepsalainen for collaboration, K. Kajantie for his 
advices and M. Laine, T. Lappi, Y. Schroder and A. Vuorinen for discussions. This work has 
been supported by the Graduate School in Particle and Nuclear Physics (GRASPANP), the 
Helsinki Institute of Physics project "Ultrarelativistic Heavy Ion Collisions" and the Academy 
of Finland, contract no. 77744. 
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Appendix A 

Pressure of the electroweak theory 



In this appendix we will write the expansion of the pressure for the full standard model. The 
detailed derivation of the result can be found in [2] and the case when the temperature is near 
the point where the electroweak crossover takes place was studied in [3]. 

To be more precise, by the full standard model we mean the theory specified by the 
Euclidean action 

£ = Jg%/?^ + ii^i^ + ^ (A.l) 
+l L Ipl L + e R Jpe R + qipqi + u R Jf)u R + d R 0d R + ig Y (q L T 2 $*t R - t R {^T 2 q L ^j , 

where we use the notation familiar from Sec. [2] and we will be keeping all the group theory 
factors for SU(iV), defined by 

T F 5 ab = TrT a T b , C F <% = [T a T b ].. 

C A 5 ab = fa™ f^, d A = 5™ , d F = 6«, 1 ' ' 

explicitly in the expressions instead of evaluating them for SU(2). Also, number of fermion 
families will be denoted by n-p and the number of fundamental scalars by ns- Mixing between 
electroweak and QCD degrees of freedom is taken into account in evaluating the pressure, but 
the pure QCD contribution to the pressure, which can be obtained from [64,65,82], is left out 
in the following computations, to be added to the total pressure in the end. The total pressure 
is then given by 

p(T) = p E (T)+p M (T)+ PQCD (T)+T*0(g 6 ), (A.3) 

where the components Pe(T) and pm(T) are given below. Note that the gluon Debye mass 
m Dgiuon> entering pqcd(T), receives electroweak corrections, 

m D, g iuon = ^D.giuonlQCD -5s t2t ^ QCD) (^C F d F n F g 2 + ^n F g' 2 + d F gl^j , (A.4) 
where gluon I QCD is the pure QCD contribution. 
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A.l Contribution from the scale 2ttT 



The scale 2irT contributes to the pressure via pe and the parameters of the effective theory. 
The effective theory is defined by 

Se = J d 3 x \^G%G% + X -F l3 F l3 + (A*) f (A*) + m 2 ^^ + 



+\{D i A%) 2 + \mlAlAl + \\A{A%Atf + ^B,) 2 + \m' 2 B,B Q 
+h 3 &$A%A a + ti 3 &$B B - ^g 3 g 3 B &A%T a <f> 



(A.5) 



where the couplings and masses are given by 



9l = 9 2 T, g'i = g' 2 T, 
A 3 = AT, A A = 0( 5 4 ), 

hs = \ 9 2 T, h' 3 = \g' 2 T, 



(A.6) 



+ 



g 2 (Pei + fee + 0{e 2 )) + (p E3 + 0(e)) + 0(g & ) 



(4tt) 2 



PexX + PE S g'i + PEYg'y + PE>g' 2 + Pe v 



rp2 



g' 2 (P'ei + /fe + 0{e 2 )) + [P' E3 + 0(e)) + 0(g 



/6\ 



;'2 



(A.7) 



(A.8) 



m§ (A) 



g 2 g' 2 A 

+ T 2 [g 2 ((3 A1 + p A2 e) + g' 2 (p Bl + fee) + A(fe + fee) + <£(fe + fee) 

<7 4 . </ 4 . g 2 g' 2 . \g 2 . Xg' 2 . X 2 . 

+ j^y 0aa + j^Pbb + ^ (3ab + ^fe + ^fe + ^ fe 

g 2 g 2 Y o g' 2 g\ R g 2 g 2 Y R Ml R , & R 



(A.9) 



and the mass counterterm for m 2 is 



T 2 



5m '- (47r) 2 eV 64 



- !V + ^ + ^V 2 " ^ - \xg' 2 + 3A 2 ) (A. 10) 
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The pressure p F is parametrized as 

Pe(T) = T A a E i + g 2 a E A + g' 2 a E B + ^a E \ + gyCtEY 

+ J^y {a Aa EAA + g' 4 a EB B + {gg') 2 a E AB + >?a E \\ + \g 2 a E A\ + Xg' 2 a EB \ 

+ gy a EYY + {ggy) 2 oL E AY + {g'gY) 2 ®EBY + ^gY a EY\ 
+ (gg s ) 2 a E As + (g'g s ) 2 a E Bs + {gYg s ) 2 a EY ^ 



+ v 2 T 2 



+ 



a Ev + ( 9 aEAv + 9> aEBy + ^ a E\u + 9y a EYu) 

; a El/u + T 4 -0(g 6 ). 



(4vr) 2 

All the parameters above are given below: 

^2 



Ot E A 



7T 

45 



1 + cIa + d F n s 



1 + d F + (2 + d F )N c 



n F 



1 

144 



5 5 
Ca^a + ^C F d F n s + -C F d F (l + A c )n F 



a EB = -ere 1 2* ,ns + 



1 



1 



1 



1 + t<*f + ( - + — d F ) N t 



n F 



a E \ 



d F (d F + 1) 



144 



(A.ll) 

(A.12) 

(A.13) 

(A.14) 

(A.15) 
(A.16) 
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36 
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23C'(-1) 11 C'(-3) 



+C F T F <i F 



23 _A_ 283 1 11 C'(-l) 11 C 7 (~3) 
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n s 



1 

128 



ci F 



1 19, A 881 1 23C'(-1) 11 C'C-3) 

— | In 1 h -7 H — — - 

e 3 4vrT 180 2 ' 3 £(-1) 6 C(-3) 



+4 



+d F 



23 , A 

— In — — 

54 4vrT 



l + ^c + $(l + ^ 



25, A 
— In 



283 2 11 C'(-l) 11 C'(— 3) 
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83 147, 8 2 C'(-l) 1 C'(-3) 
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27~~~4vrT 60 135 9 ' 27 £(-1) 27 C(-3) 
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+ 



+ 



17 



dp 



1 + ^N c + -L 1 + -± 



81 

1 + -N c 
9 



16 
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dp 
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1+3^ + -^ | + 



35 32 , 
y-yln2)n F 



4 
16 



1 + 



Ac 
9 



40, A 2 176, 8 32C'(-1) 16 C'(— 3) 
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27 4vrT ^ 27 45 9 7 27 C(-l) 27 C(-3) 
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OtEAB 



a EXX 
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16 



,1 19 A 881 1 23 C'(-l) 11 C'(— 3)' 
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d F (d F + 1) 



-n s 



, A 31 1 
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= ( 4 + 4 f^+ 4hl i^)<* + i> 
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PsY = (-^ + ^ln2 + 7 + ln^C 3F iV c (A.61) 

Pxy = (-^ln2-^ 7 -^lnA;)(d F + l)iV c (A.62) 

Pyy = \l+\^^f (A-63) 
A. 2 Contribution from the scale gT 

The contribution from the soft scale gT to the pressure is given by pm 
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+ g 3 h 3 m D B Bha + g1h 3 m 3 B Ahf ] , (A.64) 
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for which the parameters are given by 



4e 4 3 2 2m 3 J \4e 3 3 2 2m 3 / 

+ c A c FdF ( A + 1 - 3 m 2 + 5 m ^ - 1 m (A.6B) 

^ 16 b V 4e 4 3 2 2m 3 / 16 \4e 3 3 2 2m 3 / 

= **«(-s-f -T +3ba -! ta ^) (A - 67) 

7?/L4<i ~ CaCf^f^S 1 -T- - - TlnT^ 7T ln : 



8e 24 4 2m D 2 2(m 3 + m D ) 

+ Cid A (_| + 2I ln2 _^ (A.68) 

B A \f = C F (( F n s (-4 + 81ii2) + C F d F (<( F + l)ns Q +3 + 61n ^) < A - 69 > 

B BX j = ^d F (d T + l)«s (j ~ 1 4-81n2 4-61n ^-J (A.70) 

Baa/ = *(* 4- l)ns f-j - 84-41n2- 61n^-j +id F (*4- l) 2 (A.71) 

B Aha = C FrfF<iA „ s (i + ? + 21„^-4- 1 „^) (A.72) 

= ^(^! +2b 2^ +in 4) < A - 73 » 

= ft*n s (i4- S + 21„|2- + l„^) (A.74) 

*». = i*^(I + | + 2lB^ + ln^) (A.75) 

B -/ = ^A<! F -,s( i j4- 3 + 2m^ + ln£-) (A.76) 

B >>" = * dA " S (-i- 4 - 2to 2^4^- ln 2^) < A77 » 

= ^ A " 3 (-^- 4 - 2 '" 2K, + mD ) '' D ^) (AJ8) 

J*V = *«s(-i-4-21„^^-ln|i-) (A.79) 

B - = *" s (^- 4 - 2ta 2(S^- ta 2^) < A ' 80 > 
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A. 3 Pressure near the electroweak crossover 



When evaluating the pressure close to the temperatures corresponding to the electroweak 
crossover, we again need to recalculate pm(T)- The computation was carried out in [3] and 
we cite the results here. The expression for the pressure will be given by 



p(T) = PE{T)+PMi(.T)+ PM2 (T)+ PQCD {T)+T 4 0(g 



5.5\ 



(A.82) 



where Pe(T) and pqcb(T) remain the same as in the high temperature calculation. We get 
Pmi{T) from Sp in Eq. (|A.5|) by integrating over the adjoint scalars Aq and Bq. The result is 



Pmi(T) 



- - 

4vr V 3 



1 



— ( -d A ml + -mg J + j^C A d A gim z D [ --=- - - A - In 



,2 2 



1 3 

4e 4 



\ 

2m D J 



1 



+ 



(47r) 3 e 
1 

(4^)3 

,2 12, 



-C A CF5|mD + d A hlm D + ti£m' D + -C F glgf(mn + m' D ) 



dp 

T 



93 m D 



11 



7T" 



CUa -777 + -r In 2 - — + C A C F d F -- - - In 



24 6 
(mo + m' D ) ( -4 - 2 In 



+ S3#3 c Fd F - 



+ ^m D d A dF I -4-31n-^_ ] 

2m D J 



1 



^3 \ 



"ID + "T-D 



+ ti£m' D d F I -4- 3 In 



7TT-D In 
M3 



2m' T 



D 



in 

2 4 2m D /. 

(A.83) 



m n In 7 

2m D 2m' 



Theory describing just the fundamental scalar and the gauge fields is defined by 

S{, = |d 3 x S^G^ + ^Fij + iD^^D^ + ml^^ + Xs^^fY (A.84) 



where the couplings are the same as the corresponding couplings in Se, (gi, £f 3 , A 
(93 1 <?3 2 5 A3). Mass of the fundamental scalar receives corrections, 



ml 



1 



1 



1 

"2^ 



47T 



dAh^m-o 1 + In 



1*3 \ , 1 



2m D / 



+ -Ji^'v [ 1 + In 



/'3 



With parameters defined as such, we get for pm2{T): 



PM2(T) 



d F 
6n 



m 3 3 



m% 



(An) 



2 1 *■ ~I2 » ' 1 



e + 0(<? 4 ). (A.85) 



d F (d F + 1)A 3 + ( Cf?3 + 7S3 ] ( - + : > + 1 I" 



_M3_ 

2m 3 , 



.(A.86) 
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